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Summary 

The theory for predicting helicopter inflow-wake 
velocities called flat wake theory was correlated with 
several sets of experimental data. The theory was 
developed by V. E. Baskin of the USSR in 1976, and 
a computer code known as DOWN was developed 
later at Princeton University to implement the the- 
ory. The theory treats the wake geometry as rigid 
without interaction between induced velocities and 
wake structure. The wake structure is assumed to 
be a flat sheet of vorticity composed of trailing ele- 
ments whose strength depends on the azimuthal and 
radial distribution of circulation on a rotor blade. 
The code predicts the three orthogonal components 
of flow velocity in the field surrounding the rotor. 
The predictions can be utilized in rotor performance 
and helicopter real-time flight-path simulation. The 
predictive capability of the coded version of the the- 
ory was correlated with flow velocity data from sev- 
eral sources. In general, the coded version of flat 
wake theory provides vertical inflow patterns simi- 
lar to experimental patterns. Several descriptions 
of rotor blade circulation were used in the code to 
bring about a better fit of prediction of vertical inflow 
to measurement, but none were wholly satisfactory. 
To some extent, the disparities may be attributed 
to rotor- hub wake effects. However, the fundamen- 
tal wake modeling of flat wake theory may not be 
adequate for improving the correlation of code pre- 
diction with experimental data. Several assumptions 
are made in the theory that may contribute to the 
disparity between prediction and data. 

Introduction 

The prediction of inflow to a helicopter rotor and 
the wake-induced velocities arc vital to the calcula- 
tion of airloads on the rotor blades of a helicopter and 
on the other components. Inflow and wake velocities 
affect the local angles of attack and, therefore, the 
airloads that the helicopter rotor and other compo- 
nents experience. Many theories and analyses have 
been developed for such predictions of both mean and 
time- varying inflow- wake patterns. Some of these are 
described in references 1-17. Computer codes to im- 
plement these predictions have grown in complexity 
as the theories have grown in complexity. Rotor per- 
formance codes use inflow models that range from 
the assumption of uniform inflow to complex, time- 
varying, vortex filament, “free- wake” models that 
compute complex inflow distributions at the rotor 
disk. These generally ignore the effects of the fuse- 
lage on the velocities seen by elements of the lifting 
rotor. Even though computer processing time and as- 
sociated costs have dramatically decreased in recent 


years, there is still a need for codes that are simple 
and fast and , at the same time, offer accurate repre- 
sentations of inflow" and wake patterns (as mentioned 
in refs. 12, 17, and 18). 

“Real-time” simulation of helicopter flight for 
agility and controllability studies requires faster and 
more representative inflow- w r ake modeling (ref. 17). 
Reference 18 concludes that an understanding of the 
influence of the rotor wake on the fixed tail surfaces 
and tail rotor is one of the most important factors 
needed to improve the representation of helicopter 
flight path motions. Reference 19, which describes 
the development history of the U.S. Army’s AH-64 
helicopter, provides examples of the importance of 
the w^ake interaction between the main rotor, tail ro- 
tor, and empennage. The empennage design had to 
be changed as a better understanding of wake effects 
evolved. 

The calculation of dynamic rotor characteristics 
has depended generally on the assumption of uni- 
form inflow (invariant with azimuth or radial posi- 
tion). This, of course, is unrealistic. As early as 
1934, Wheatley concluded in reference 1 that “The 
blade motion is critically dependent upon the distri- 
bution of induced velocities over the rotor disk and 
cannot be calculated rigorously without the accurate 
determination of the induced flow.” In recent years 
more sophisticated (and complex) inflow descriptions 
have come into use, but as noted in reference 17, 
these have not generally resulted in significantly bet- 
ter predictions. 

A theory that may offer an acceptable balance 
of speed and prediction capability for mean inflow- 
wake velocity has been used in Eastern Europe and 
is referred to as the flat wake theory (ref. 20). The 
theory and associated assumptions and limitations 
are well-described in references 20, 21, and 22. The 
theory involves assuming a rigid undeflccted ribbon 
of vorticit}' generated by the rotor blades as shown 
in figure 1. The vortices transferred from the ro- 
tor blades to the slipstream of the rotor are as- 
sumed to merge to a single ribbon of vorticity trail- 
ing behind the rotor. At advance ratios greater than 
0.13 (52 knots for a typical rotor system), the non- 
linearities of the vortex motion arc judged in refer- 
ence 20 to have a minor influence on the induced 
velocities near the rotor disk plane. The flat w r akc 
theory is most useful for predicting time-averaged in- 
flow and wake characteristics about a rotor. Since the 
prediction is only from a time- aver aged sense, it is as- 
sumed to be unaffected by blade passage. For a time- 
varying rotor flow field analysis, references 11-15 can 
provide an estimate of the influence of blade passage. 



Such theories are useful for time-varying blade loads, 
rotor dynamics analyses, and rotor acoustics. 

The results in reference 20 indicate good agree- 
ment with some experimental measurements of the 
velocity components in the wake of a lifting rotor. 
The details of the calculation procedure developed 
for reference 20 were unavailable, and the data used 
for validation were limited to wake velocity measure- 
ments (not inflow); therefore, a study of the flat wake 
theory and its possible application as a computer 
code with an acceptable balance of speed and pre- 
diction capability w r as initiated. Using the guidelines 
offered in references 20 and 21, a relatively compact 
computer code known as DOWN for the implemen- 
tation of the flat wake theory was developed in recent 
years by Howard C. Curtiss and Robert M. McKillup 
of Princeton University with a NASA research grant. 
A user’s guide (ref. 23) is available for the code. 

Predictions of inflow and wake velocities from 
DOWN are compared with the data of several 
sources. Until recent years there has been very lim- 
ited experimental data describing rotor inflow and 
wake (mostly in ref. 5). This deficiency has been 
partially remedied with the data presented in refer- 
ences 24 -34. These data consist of measurements 
of the longitudinal and vertical velocity increments 
caused by a rotor for a variety of advance ratios, rotor 
blade planforms, and measurement planes. The data 
of references 24-34 have been used for the evaluation 
of other inflow- wake analyses (refs. 11, 35, and 36). 
Some of these data, along with what is available de- 
scribing rotor flow field characteristics in reference 5, 
offer a means for the calibration of DOWN. The ver- 
tical induced velocity component is of greatest inter- 
est since it is critical in the calculation procedures for 
rotor performance analyses. The lateral component 
of induced velocity affects sidewash flow angles in 
the region of a helicopter vertical fin and tail rotor, 
thus affecting helicopter yaw stability and control. 
The limited available data in reference 37 describing 
lateral velocities are also used in the calibration of 
DOWN. 

Symbols 

The coordinate system is shown in figure 2. Note 
that this convention differs from the one used in the 
United States. 

Cp rotor thrust coefficient, 

T/[p a (nR ) 2 nR 2 } 

Gl, G2, G3 assumed circulation dis- 

tribution on rotor blade 
without azimuthal variation 
(correction), T f (p) 


Glm, G2m, G3m 

assumed circulation distri- 
bution on rotor blade with 
azimuthal variation (correc- 
tion), r'(p) + r'{p,ip) 

R 

rotor radius, ft 

r 

rotor blade station, ft 

A r 

increment in rotor blade 
station, ft 

T 

rotor thrust, lb 

V 

flight speed, ft/sec or knots 

v t 

rotor tip speed, fif?, ft /sec 

V x 

increment in longitudinal 


velocity (positive forward) 



due to rotor inflow and 
wake, ft /sec 

v y 

increment in vertical veloc- 
ity (positive up) due to ro- 
tor inflow and wake, ft /sec 

V z 

increment in sidewash 
velocity (positive right) due 
to rotor inflow- and wake, 
ft /sec 

x-.y,z 

Cartesian coordinate system 
centered in rotor hub 
with x positive forward, 
y positive upw-ard, and z 
positive laterally to right 
(convention used in ref. 20 
whereby xz is the tip path 
plane); nondimensionalized 
by radius, R (see fig. 2) 

a 

angle of attack of rotor, deg 

r 

relative circulation averaged 
around rotor, ft 2 /sec 

r' 

blade circulation, ft 2 /sec 

Ar 

increment in blade circu- 
lation as a function of az- 
imuth, ft 2 /sec 

Ar, 

increment in sidewush 
velocity ratio, v z /Vt 

AX 

increment in rotor inflow 
ratio, VyjVt 

R 

advance ratio, V/Vt 

P 

rotor blade station, r/R 

A p 

increment in rotor radius, 
Ar/R 
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Pa air density, slugs/ft 3 

ip blade azimuth position, deg 

£7 rotor rotational speed, 

radians/sec 

Discussion of Theory and Computer 
Code 

Theory 

The development of the flat wake theory and an 
outline of the means to implement it are described 
in references 20 and 21. It begins with the consid- 
eration of an idealized rotor vortex system that is 
formed by shed and trailing vortices formed by a 
rotor blade in hover or axial flight. These form a 
sheet of vorticity of a classical helical shape called 
the vortex sheet behind each blade (fig. 3). In ax- 
ial flight, i.e., hover or vertical climb, the sheets of 
several rotor blades appear to form a column. The 
column has an axis that is perpendicular to the rotor 
disk. In horizontal flight an inplane component of 
velocity results and the axis of the column is skewed 
with respect to the rotor. As the horizontal flight 
velocity component increases, the column of vortic- 
ity below the rotor deflects more to the rear. At 
sufficiently high flight speeds, the column becomes 
practically flat and results in a “ribbon” of vorticity 
transported downstream at the free-stream velocity 
(fig. 1). Within this ribbon the trailing vortices leav- 
ing the blade have a cycloidal form (fig. 4(a)). These 
trailing vortex elements are chosen to start at several 
locations along the blade span and result in a trajec- 
tory that suggests a weaving of the vortices within 
the ribbon. The ribbon is assumed to be a “rigid” 
wake, and therefore there is no deformation of the 
wake structure because of wake induced velocities. 

The trailing vortices that are formed are the 
consequence of radial and azimuthal distributions 
of circulation on the rotor blade. The variation of 
circulation along the blade span leads to the presence 
of trailing vortices springing from the blade. As a 
first approximation, Baskin (ref. 20) assumes that 
the variation of circulation with azimuth may be 
neglected for practical calculations if, at each blade 
station, the circulation is averaged over a complete 
rotor revolution where in 

r(/>, /x, i/>) « r'(p) + at (n, ip) 

we have 

AT(p,ip) = 0 

However, he does offer a generalization to this first 
approximation to account for azimuthal variation. 


Neglecting higher order terms, this becomes 

AT(/^, ip) ~ — 1.5/x sin ip + 1.12/i 2 (l — cos 2 ip) 

He also assumes that the bound vorticity contribu- 
tion to velocity increments can be ignored. He states 
that bound vortices having the same circulation at 
all azimuth angles would not generate any vertical 
induced velocities. However, they would produce ro- 
tational induced velocities at points located above as 
well as below the rotor disk. On the average, the 
rotational induced (inplane) velocity represents no 
more than 0.5 percent of the blade tip speed and can 
be neglected. 

Although this classical helical wake representa- 
tion has been shown to be a good representation of 
the physical characteristics of the rotor wake, math- 
ematical implementation of this geometry can be 
quite extensive and too complex for predesign and 
real-time simulation analytical methods. Therefore, 
Baskin has developed an equivalent rectilinear wake 
structure involving lateral and longitudinal vortices. 
(See fig. 4.) The rotor system is represented as a 
lifting disk with a lattice of vortices distributed in a 
rectilinear fashion. Radial and azimuthal variations 
in rotor circulation over the disk and their effect in 
the shed wake can be accounted for at each of the 
“nodes” of the lattice. The wake is represented by a 
system of longitudinal trailing vortices from each of 
the lateral lattice locations (fig. 4(b)). No transverse 
or shed vorticity is modeled in the wake. 

In computing the resultant induced velocity, the 
velocity components generated by various vortex sub- 
systems forming the wake can be computed sepa- 
rately and then superimposed. The Biot-Savart law 
is used in the calculation of the field-point velocities 
induced from each of the lateral and longitudinal vor- 
tex segments. Reference 21 offers a detailed descrip- 
tion of the contributions of the various vortices to 
the v x , v y , and v z components and of the process of 
summation. 

Experience and theoretical considerations cited in 
reference 20 indicate that this flat wake of vortic- 
ity is applicable for an advance ratio p greater than 
1.63 yfUj" For a typical rotor tip speed of 700 ft/sec 
and Cx = 0.006, the corresponding advance ratio is 
approximately 0.13 and the flight speed is 52 knots. 
In reference 16 Ormiston agrees that “A simple, flat 
planar wake has been found to be a valid configura- 
tion for the wake vorticity down to advance ratios of 
approximately 0.15 for a nominal rotor thrust coef- 
ficient.” If a sinusoidal variation of rotor blade cir- 
culation with azimuth angle is allowed (as suggested 
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by refs. 6, 16, and 20), the flat wake analysis may be 
applicable at a moderately lower advance ratio. 

Computer Code 

The study and computer coding of flat wake the- 
ory of reference 20 was arranged with a NASA re- 
search grant to Princeton University. The computer 
code known as DOWN that evolved can be used to 
calculate mean induced flow field velocities (vertical, 
longitudinal, and lateral) with few input parameters 
(/x, Cti and location of prediction plane) for rotor 
operating conditions. A selected circulation distribu- 
tion on the rotor blade as a function of radius must 
be defined that is to represent an average for 1 revo- 
lution. These are identified in the figures by Gl, G2, 
or G3. Here, Gl represents a linear distribution, G2 
represents a modified “parabolic” distribution, and 
G3 represents a modified “quartic” distribution. The 
respective equations are as follows: 

For Gl, 

T' = 10.0(0.0589 + 1.3783/?) 

for G2, 

r' = io.o(p 2 - p 3 ) 

and for G3, 

r' = [(1.0 - 0.685 p)p(p 2 + 0.011)]/(p 2 + 0.0055) 

Figure 5 shows these circulation distributions as a 
function of blade radius where the patterns are car- 
ried into the hub and the blade root cutout is ig- 
nored. The quartic distribution is more representa- 
tive of those shown in reference 38 than of either 
the linear or parabolic distributions. The user of the 
program DOWN can easily introduce alternate dis- 
tributions and also has the option of implementing 
the azimuthal correction suggested by Baskin in ref- 
erence 20. When the azimuthal correction is used in 
this paper, the circulation distributions are identified 
as Glm, G2m, and G3m. 

The circulation distribution on a rotor blade is 
normalized within the code as follows: 

, r'M 

2 firmed, 

so that any distribution may be used. There is then 
a summation of incremental velocity components (in 
the Cartesian coordinate system) that are a function 
of the incremental change of circulation distribution 
spanwise (in the polar coordinate system) on the 
rotor blade. Within the code the increment in radius 
can be changed. Generally, a selection of increment 


in radius A p of 0.1 and the associated increment 
(or decrement) in circulation strength are adequate. 
Using smaller increments does improve accuracy, but 
experience cited in reference 20 indicates that the 
gain is not significant. Both the “reversed flow 
region” and the local blade stalled angle of attack 
due to interaction with blade tip vortices are not 
accounted for; also the influence of the wake of the 
rotor hub is not accounted for. However, these effects 
are typically ignored in other simple inflow analyses 
as well in more complex methods (ref. 17). 

The code provides tabulations of induced veloc- 
ity components normalized with tip speed in either 
Cartesian or cylindrical coordinate systems with two 
dependent coordinates. The axes convention of ref- 
erence 20 is maintained (fig. 2) even though it differs 
from the U.S. convention. 

Experimental Data 

The data that are used for the calibration of 
the flat wake code are described in the order used 
for correlation, and the test conditions are listed in 
table I. 

Experimental Data From Reference 20 

The data presented in reference 20 are mean 
wake velocity measurements for an advance ratio of 
0.15 and a thrust coefficient of 0.006. There is no 
other information regarding rotor angle of attack, 
power, propulsive force, or data accuracy. The test 
apparatus and flow measurement technique used are 
not described. The rotor appears to have had tapered 
planform blades. The data, which are mean vertical 
wake velocities nondimensionalized with rotor tip 
speed, are rather limited for an evaluation of the 
original version of the flat wake analysis. Baskin 
was able to obtain a good fit of flat-wake-theory 
prediction to these data, as shown in figure 6. 

Experimental Data From 

References 24-34 

The data presented in references 24-34 are laser 
velocimetry measurements of longitudinal and verti- 
cal velocities one blade chord or less above the ro- 
tor tip path plane. Test conditions of advance ratio, 
thrust coefficient, rotor angle of attack, and measure- 
ment plane height are listed in table I. Other rotor 
operating conditions such as power, propulsive force, 
etc., can be found in the references. Velocity mea- 
surement accuracy was defined in references 24 34 as 
approximately 2 percent, which translates into an ac- 
curacy for AA ranging from approximately ±0.0002 
at fi ~ 0.149 to ±0.0008 at ji — 0.400. The measure- 
ment locations were made at radial stations from 
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p — 0.2 to p — 1.2 at every 30° of azimuth. The 
measurements are nondimensionalized with rotor tip 
speed Vt for both mean and time-varying velocities. 
For the purpose of evaluating the flat wake code, only 
the mean data from references 24, 25, 26, 29, and 30 
were utilized for this report. These provide the data 
obtained in a plane parallel to and 1 chord above the 
rotor tip path plane for rectangular planform blades 
operated at advance ratios of 0.15, 0.23, 0.30, 0.35, 
and 0.40, respectively. 

These data were obtained with a two-component 
laser velocimeter (LV) system operated in the NASA 
Langley 14- by 22-Foot Subsonic Tunnel (ref. 39). 
The LV system measured the instantaneous compo- 
nents of velocity in the longitudinal and vertical di- 
rections. A general-purpose rotorcraft model system 
called the 2-meter rotor test system (2MRTS) was 
part of the experimental apparatus. A more detailed 
description of the 2MRTS is provided in reference 40. 
The rotor model provided a reasonable representa- 
tion of rotor, hub, and fuselage so that the hub and 
the presence of the fuselage affect flow velocities in a 
manner similar to that on full-scale helicopters. 

Experimental Data From Reference 5 

In reference 5, Hey son and Katzoff described a 
theoretical and experimental investigation of the ve- 
locities near a lifting rotor. The analysis and data 
have been used for many years for helicopter de- 
sign. Mean flow velocity measurements were made 
above and below the rotor using five-hole pressure- 
measuring probes. Velocity measurement accuracy 
was 4 percent with a corresponding accuracy of AA 
of ±0.0015 for the data used herein. The model had a 
teetering rotor system with two untapered, untwisted 
blades with NACA 0012 airfoils. The rotor test appa- 
ratus had no fuselage where the rotor was driven by 
a shaft enclosed within a tall cylindrical windshield. 
Thrust coefficients ranged from 0.00321 to 0.00373, 
tip speed ranged from 450 to 500 ft/sec, and disk 
loading ranged from 1.55 to 2.21 lb/sq ft. Data were 
obtained for advance ratios of 0.095, 0.139, 0.140, 
and 0.232. 

Thrust coefficients, tip speed, and disk loading 
were quite low compared with typical conditions 
for current helicopters. Both test and analysis of 


reference 5 deal with lightly loaded rotors and, as 
pointed out in reference 3, may pose serious limita- 
tions to the results. As a cylinder, the windshield 
may have had a greater effect on flow velocities than 
would be the case for a fuselage representation as in 
references 24-34. It is acknowledged in reference 5 
that the windshield affected the flow patterns. Al- 
though most data and analysis of reference 5 were for 
rotor and test conditions that are not representative 
of current rotors and flight conditions, the informa- 
tion for p = 0.232 is useful for the calibration of the 
flat- wake- theory code. 

Experimental Data From Reference 37 

The primary purpose of the wind tunnel program 
described in reference 37 was to gain an understand- 
ing of the flow field environment just ahead of the 
empennage system as affected by the hub and ro- 
tor wake. The fuselage was a model of a U.S. Army 
attack helicopter. Wake surveys were made that de- 
fined the downwash and side wash angles and the dy- 
namic pressure ratios in that region. The surveys 
were made at four vertical heights, but since the rotor 
used in the test was undersized, the effective longi- 
tudinal location was x = —1.76. (A typical helicopter 
vertical tail position would be x w —1.20.) The model 
fuselage was 30 percent scale, whereas the rotor was 
21 percent scale. Wake measurements were obtained 
both with and without rotor blades so that hub wake 
effects could be defined. Since the hub was not re- 
moved, a definitive appraisal of purely hub effects 
could not be made. 

Of the test conditions, only one was useful for 
the calibration of the flat-wake-code prediction of 
sidewash velocity, and that was at p = 0.29. The rest 
were at very low advance ratios. It was noted that 
the rotor wake was strongly asymmetric where the 
downwash behind the advancing side of the rotor 
was greater than behind the retreating side. The 
rotor conditions were Cj = 0.006 and a = —4.2.° 
The rotor propulsive force and torque coefficients 
arc not provided in the reference, but the rotor was 
trimmed to zero flapping. The flow survey made with 
only the hub operating (i.e., no blades) showed that 
the hub contributes a substantial component of the 
lateral velocity in the vicinity of the vertical tail. 
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Presentation of Results 


The results of this evaluation are presented in the following figures: 

Figure 


Circulation distributions used in predictions . 5 

Data and predictions from reference 20; AA versus ip 6 

DOWN predictions compared with reference 20 data: 

Data from reference 20 with Gl, Glm, G2, G2m, G3, and G3m; AA versus ip . . 7 

Data from reference 20 with G2 and G3; A A versus p 8 

Data from reference 20 with G3 and G3m; A A versus 9 

DOWN predictions of AA versus ip for p — 0.60, 0.79, and 0.98: 

With data from reference 24; p = 0.149 10 

With data from reference 25; p — 0.230 11 

With data from reference 26; p = 0.300 12 

With data from reference 29; p — 0.350 .13 

With data from reference 30; p = 0.400 14 

Fuselage- corrected A A for references 24, 25, and 26; A A versus ip 15 

Application of fuselage-corrected A A to reference 25 data; A A versus ip . . 16 

“Three-dimensional” plots of DOWN predictions and data; A A versus x and z: 

Data from reference 24; p = 0-149 17 

Data from reference 25; p = 0.230 18 

Data from reference 26; p = 0.300 19 

Data from reference 29; p = 0.350 20 

Data from reference 30; p = 0.400 21 

DOWN predictions compared with data from reference 5; AA versus z: 

At x = 1.0, 0.5,0, -0.50, and — 1.07; y — 0.16 to 0.19 22 

At y = 0.19, —0.07, and —0.21; x = —0.50 23 

At y = 0.21, -0.06, and — 0.19;x = —1.07 24 

DOWN predictions compared with data from reference 37: 

At? versus z . . 25 

AA versus z 26 


Discussion of Results 

Correlation With Data From Reference 20 

The calculated induced velocities of program 
DOWN should correspond to the analytical predic- 
tions of reference 20 since DOWN was meant to be 
the implementation of the analysis of reference 20. 
Figure 6 shows both data and predictions reproduced 
from reference 20. Figure 6 shows the azimuthal dis- 
tribution of induced vertical velocity A A at p = 0.70, 
p = 0.15, Ct = 0.006, and y = —0.10. According to 
reference 20, figure 6 shows satisfactory agreement 
between prediction and experimental data. It is also 
judged in reference 20 that using a circulation distri- 
bution that does not account for azimuthal variation 


fits better to the data than a circulation distribution 
that is a function _pf radius and azimuth.. 

Figure 7 shows a comparison of DOWN pre- 
dictions (utilizing the three different circulations, 
with and without the azimuthal variation) with ref- 
erence 20 data. DOWN predictions are for the 
linear circulation (fig. 7(a)), parabolic circulation 
(fig. 7(b)), and quartic circulation (fig. 7(c)). Of 
the distributions, the quartic distribution seems to 
provide the best fit, though only a fair fit, whereas 
the linear distribution fits best in the azimuth range 
from 135° to 240°. In the azimuth range from 330° to 
30°, all predictions and those of reference 20 (shown 
in fig. 6) show peaks that do not fit well with the 
data. They likely show the effects of model mounting 
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apparatus on flow velocities. The azimuthal variation 
of circulation distribution does not improve the fit to 
the data. 

Figures 8 and 9 present data from reference 20 
where test conditions were p = 0.15, C? = 0.006, 
and y — —0.10. The induced velocity AA is plotted 
against p for various azimuth combinations. Where 
azimuths are greater than 180°, a negative range 
of p is used. The azimuth combinations of 0° and 
180° (figs. 8(a) and 9(a)), 45° and 225° (figs. 8(b) 
and 9(b)), 90° and 270° (figs. 8(c) and 9(c)), and 135° 
and 315° (figs. 8(d) and 9(d)) are shown. Figure 8 
compares DOWN predictions having the parabolic 
and quartic circulations with data from reference 20. 
Figure 9 compares DOWN predictions having the 
quartic circulation (with and without the azimuthal 
variation) with data from reference 20. (Recall that 
the notation m in the figures (e.g., G3m) indicates 
that an azimuthal variation of circulation is used in 
the prediction.) The results of DOWN for parabolic 
and quartic circulations appear to fit to the data, 
with the parabolic circulation being somewhat better 
in figure 8 than the quartic circulation. However, the 
fit is no better than that of figure 7. In figure 9, 
the use of an azimuthal variation of the quartic 
circulation did not improve the fit of prediction to 
data. 

Figures 7, 8, and 9 do not show a clearly bet- 
ter correlation of DOWN predictions with one or the 
other of the circulation descriptions from the data 
of reference 20. This analysis was unable to repro- 
duce the “good” agreement obtained in reference 20 
(shown in fig. 6). Since the coding of DOWN is 
based on the guidance of reference 20, a closer fit 
of the DOWN prediction to the analysis of refer- 
ence 20 should be expected. The results shown in 
reference 20 cannot be verified, whereas the code 
DOWN has been carefully checked. The poorest fit 
is in the azimuth range where rotor-hub wake effects 
must be significant. 

Correlation With Data From 

References 24, 25, 26, 29, and 30 

Comparisons of vertical induced velocity data 
of references 24, 25, 26, 29, and 30 (volumes I, 
II, III, VI, and VII, respectively) were made with 
the DOWN predictions and are presented in fig- 
ures 10-14. Shown is the variation of induced flow 
AA with azimuth xp for three selected radial stations 
(p = 0.60,0.79, and 0.98). The data shown have 
been corrected for the influence of the fuselage as 
defined in references 41 and 42. This correction 
process is demonstrated in figures 15 and 16 for 
p = 0.15,0.23, and 0.30 and at radial stations of 


0.60, 0.79, and 0.98. Figure 16 suggests that the 
experimental measurements are modestly affected by 
the presence of the fuselage as suggested by the 
analyses of reference 41. 

Since there was no clear “best” circulation de- 
scription for correlation with experimental data in 
figures 7, 8, and 9, and the quartic distribution 
appeared to be similar to patterns shown in refer- 
ence 38, the quartic distribution was chosen for cor- 
relation of DOWN predictions with the data from 
references 24, 25, 26, 29, and 30. In general, the 
code provides a prediction of vertical inflow charac- 
teristics similar to that of experimental data as shown 
in figures 10-14. The data for the radial station 0.98 
fitted best for ^ = 90° to 270°, but with a signif- 
icant difference in the azimuth range from 315° to 
60°. That characteristic appears for the other ra- 
dial stations of 0.60 and 0.79 as well (and in ref. 20). 
There are some other differences between predictions 
and data for these two stations. For example, in fig- 
ures 11(a), 11(b), 12(a), and 12(b) there is signifi- 
cant upwash shown for the azimuth range from 60° to 
150°, whereas the prediction shows downwash. The 
inaccuracy of velocity measurement is not enough to 
explain the difference. The hub wake effect on flow 
patterns to the rear and below the rotor would be 
even more significant. There are no data obtained 
in the tests cited in references 24 to 34 to provide 
an appraisal of such a hub wake effect. Even re- 
cently developed analytical tools (such as those of 
refs. 43 and 44) for rotorcraft interactional aerody- 
namics have yet to address that need. It does ap- 
pear that accounting for azimuthal variation of the 
quartic circulation offers a modest improvement in 
correlation. 

Figures 17-21 offer another way of comparing 
DOWN predictions with the data from references 24, 
25, 26, 29, and 30, i.e., a three-dimensional format. 
The differences between prediction and data can pro- 
vide an appreciation of the complete inflow as op- 
posed to those shown in figures 10-14. An obvious 
difference is the upwash predicted on both sides of 
the rotor disk for a radius where p > 1.0 and for az- 
imuth ranges from 60° to 90° and 270° to 300° which 
are not evident in the data. However, the predicted 
patterns do correlate with data in showing maximum 
downwash occurring in the rear portion of the rotor 
disk and skewed to the advancing blade side of the 
disk (azimuth range of 0° to 90°). That was not the 
case in a comparison of the predictions of the analyt- 
ical method of reference 14 with the very same data, 
as concluded in reference 35. The results of refer- 
ence 35 using CAMRAD (Comprehensive Analytical 
Model of Rotorcraft Aerodynamics and Dynamics) 
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have shown maximum downwash toward the retreat- 
ing side of the rotor disk (ref. 14). Several other 
analytical methods reviewed in reference 36 had the 
same problem of predicting the largest values of in- 
duced inflow to be on the retreating side of the rotor 
disk, which is the opposite from that measured. So, 
it appears that DOWN is better than those analyses 
(CAMRAD, UTRC Rotorcraft Wake Analysis, and 
Beddoes) mentioned in references 35 and 36 in pre- 
dicting the largest values of inflow on the retreating 
side of the rotor disk. 

Correlation With Data From Reference 5 

Several sets of data (all at fi = 0.232) presented 
in reference 5 are used for correlation. The first set 
(fig. 22) provides for a correlation with the prediction 
of inflow in a plane y « 0.18 above the rotor tip path 
plane for ji = 0.232. The several longitudinal stations 
shown are at x = 1.0, 0.5, 0, —0.5, and —1.07. The 
measurement height at several stations differed mod- 
estly from 0.18, but the correct height is noted on 
each plot. The predicted inflow patterns are in gen- 
eral agreement with the data although there is a 
greater downflow measured for z < —0.5. Figures 23 
and 24 offer comparisons of inflow and wake char- 
acteristics, i.e., above and below the rotor at two 
longitudinal stations (x = —0.50 and —1.07, respec- 
tively). The A A characteristics shown in figure 24 
are at a longitudinal station representative of an em- 
pennage station (x — —1.07). Again, the predicted 
airflow characteristics are in general agreement with 
the data, but there is a greater disparity for 2 < 0. 

Correlation With Data From Reference 37 

The wake survey test data of reference 37 may be 
more representative of the airflow characteristic of 
a helicopter empennage location than that of refer- 
ence 5 since there is a full representation of a fuselage 
with the resultant wake. The results of reference 37 
indicate that the flow field is quite chaotic with a 
strong swirling flow, which appears to be mainly the 
influence of the rotor hub. The hub- alone data were 
subtracted from data for the rotor to obtain wake 
flow characteristics more appropriate for correlation 
(since flat wake theory does not address hub effects). 
In figures 25 (for sidewash A rj) and 26 (for down- 
wash A A) the rotors with and without hub data are 
compared with the flat-wakc-code predictions. The 
agreement between data and prediction is very poor. 
A reason may be that separated flow from the fuse- 
lage and engine nacelles is substantial and the resul- 
tant wake has a gross effect on patterns at the tail 
pylon. It is obvious that just subtracting hub-alone 
data from rotor-plus-hub data does not offer much 


improvement in correlation and there is a need for 
more sophisticated analyses, such as the free wake 
flow field analysis of reference 44 wherein a better 
correlation with data from reference 37 is achieved. 
However, the good correlation in reference 44 can- 
not be attributed to the incorporation of a hub and 
fuselage wake. 

Concluding Remarks I 

An analysis of helicopter inflow-wake velocities I 

called flat wake theory was implemented in a com- j 

puter code. This code known as DOWN will pre- \ 

diet the three orthogonal components of flow veloc- I 

ity in the field surrounding the rotor. The predic- | 

tive capability of the coded versions of the theory - 

was correlated with flow velocity data from several : 

sources. These provide mostly vertical induced flow 
velocity measurements and some meager lateral in- : 
duced flow information. The data used from these s 
several sources are compromised by the significant \ 

influence of the hub and rotor drive system. 1 

The assumptions that were made in the develop- s 

ment of the flat rotor flow field analysis may well : 

contribute to a disparity between data and predic- ~ 

tion, and these are listed as follows: (1) the model- \ 

ing of the classical helical vortex wake by a system of j 

rectilinear vortices, (2) a rigid wake without interac- l 

tion between induced velocities and wake structure, : 

(3) an azimuthal average of circulation distribution 
on a rotor blade, (4) an azimuthal variation of circu- 
lation (and therefore no shed vorticity) not directly I 
incorporated in the analysis, and (5) “bound” vortic- ! 
ity effects not represented. | 

In general, the coded version of flat wake the- 
ory provides vertical inflow patterns similar to ex- I 
perimental time-averaged measurements. Since the ! 

computer code was meant to be an implementation 
of an existing calculation scheme, it at least should 
have duplicated rotor flow calculations made for the 
calibration of the original presentation of the theory. 

That was not accomplished with the code DOWN. 

Since the code has been thoroughly checked by the 
developers, the original presentation of calibration 
needs verification. 

Although several rotor blade circulation descrip- 
tions were used in the code in an attempt to bring I 

about a better fit of prediction to measured verti- ! 

cal flow velocities, none were wholly satisfactory. Of 
these circulation patterns the one that is most like a 
flight test distribution is marginally best. Applying 
a factorto represent an azimuthal variation of circu- 
lation provided a modest improvement to the fit of 
predictions to some data. To some extent the dispar- 
ities between prediction and data may be attributed 
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to rotor hub and fuselage wake effects that are not 
accounted for in the theory. 

For the prediction of induced sidewash velocities 
in the vicinity of the helicopter tail, the paucity 
of data limited the attempted correlation of the 
code. The code predictions did not fit well with the 
data. However, the flow in that region is obviously 
dominated by the swirl effected by the rotor hub and 
the separated flow (i.e., wake) from the hub, pylon, 
and fuselage. 

fn spite of the evident disparities between data 
and DOWN predictions, the code predictions may 
be no worse than those of other analytical codes. 
For example, DOWN does predict greater downwash 
velocities on the advancing side of the rotor disk than 
on the retreating side which is not done by other 
analytical codes. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
March 5, 1992 
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Table I. Test Conditions of Experimental Data 


(a) Data from references 20, 24, 25, 26, 29, and 30 


Reference 

M 

Cj' 

y 

P 

ip, deg 

a, deg 

20 

0.150 

0.0060 

-0.100 

0.7 

Range 

(a) 

24 

.149 

.0063 

.077 

Range 



-3.00 

25 

.230 

.0065 







-3.04 

26 

.300 

.0065 







-4.04 

29 

.350 

.0064 







-5.70 

30 

.400 

.0064 

' 

1* 




' 

-5.70 


a Not applicable. 


(b) Data from references 5 and 37 


Reference 


C T 

y 

X 

z 

a, deg 

5 

0.232 

0.00321 

0.19/— 0.23 

1.0, 0.5, 0, 
-0.5, -1.07 

Range 

-9.50 

37 

0.290 

0.006 

—0.22/— 0.51 

-1.76 

Range 

-4.20 





1 





Figure 1. Flat wake formed by tip vortices of 
a rotor. 



Figure 2. Coordinate system maintained from 
reference 20. Note that it differs from 
convention used in the U.S. 
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Blade radius, p 

Figure 5. Linear (Gl), parabolic (G2), and quartic (G3) circulation distributions used in DOWN for flow field 
predictions. 



yr, deg 


Figure 6. Data and flat wake predictions taken from figure 2.31 of reference 20 with circulation as a function 
of p, r'(p), and as a function of p and ip, T (p,p). p — 0.15; C T = 0.006; y = -0.10; p = 0.70. 




0 30 60 90 120 150 180 210 240 270 300 330 360 

V, deg 

(a) Linear circulation with and without azimuthal variation. 



V, deg 

(b) Parabolic circulation with and without azimuthal variation. 

Figure 7. Correlation of DOWN prediction with various circulation distributions of AA with data of refer 
ence 20. u = 0.15; Cj 1 — 0.006; p = 0.7; y — —0.10. 






(a) xp = 0° and 180°. 



P 


(b) ip = 45° and 225°. 

Figure 8. Correlation of DOWN prediction of AA with data of reference 20 by using the parabolic and quartic 
circulations for various azimuthal angles. // = 0.15; Ct — 0.006; y — —0.10. 
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(a) p = 0.60. 



(b) p = 0.79. 

Figure 10. Correlation of DOWN prediction of AA with data of reference 24 by using quartic circulation with 
and without azimuthal variation of circulation, p = 0.149; Cx — 0.0063; y = 0.077. 










W deg 

(a) p = 0.60. 



Figure 11. Correlation of DOWN prediction of A A with data of reference 25 by using the quartic circulation 
distribution with and without azimuthal variation of circulation, p — 0.230; Cj — 0.0065; y = 0.077. 
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(a) p = 0.60. 



V, deg 

(b) p = 0.79. 


Figure 12. Correlation of DOWN prediction of AA with data of reference 26 by using the quartic circulation 
distribution with and without azimuthal variation of circulation, p = 0.300; Cx — 0.00649; y — 0.077. 
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(b) p = 0.79. 


Figure 13. Correlation of DOWN prediction of AA with data of reference 29 by using the quartic circulation 
distribution with and without azimuthal variation of circulation, p = 0.350; Cj = 0.0064; y = 0.077. 














(c) p = 0.98. 
Figure 14. Concluded. 
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(b) p = 0.79. 


Figure 16. Effect of applying a correction for the presence of the fuselage to AA data of reference 25 









(a) DOWN prediction with G3m. 



(b) Data from reference 24. 

Figure 17. Comparison of DOWN prediction of AA with data from reference 24 with ^ = 0.149. 
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Figure 18. Comparison of DOWN prediction of AA with data from reference 25 with pi = 0.230. 
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(b) Data from reference 26. 
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Figure 19. Comparison of DOWN prediction of AA with data from reference 26 with // 













Figure 21. 
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(a) DOWN prediction with G3m. 
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Comparison of DOWN prediction of AA with data from reference 30 with [i 


= 0.400. 
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Figure 23. Correlation of DOWN prediction (using quart ic circulation distribution with and without azimuthal 
variation of circulation) of AA with data from reference 5 for several vertical stations at x = —0.50. 
u = 0.232; C T = 0.00321. 
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Figure 24. Correlation of DOWN prediction (using quartic circulation distribution with and without azimuthal 
variation of circulation) of AA with data from reference 5 for several vertical stations at x = —1.07 
H = 0.232; C T = 0.00321. 
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(a) y = -0.21. 
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(b) y = -0.31. 

Figure 25. Correlation of DOWN prediction (using quartic circulation distribution) of At? with data from 
reference 37 for several vertical stations at x = -1.70. y = 0.29; Cj = 0.006. 
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(a) y = -0.21. 



(b) y = -0.31. 

Figure 26. Correlation of DOWN prediction (using quartic circulation) of AA with data from reference 37 for 
several vertical stations at x — —1.70. n — 0.29; Cj — 0.006. 
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